Tunneling Spectroscopy of Two-level Systems Inside Josephson Junction 
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We consider a two-level (TL) system with energy level separation h£lo inside a Josephson junction. 
The junction is shunted by a resistor R and is current I (or voltage V = RL) biased. If the TL 
system modulates the Josephson energy and/or is optically active, it is Rabi driven by the Josephson 
oscillations in the running phase regime near the resonance 2eV = hQo- The Rabi oscillations, in 
turn, translate into oscillations of current and voltage which can be detected in noise measurements. 
This effect provides an option to fully characterize the TL systems and to find the TL's contribution 
to the decoherence when the junction is used as a qubit. 
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Josephson junctions are promising candidates for 
qubits in quantum computing. However, to ensure long 
decoherence time, the phase difference across the junc- 
tion should be coupled very weakly with other low- 
energy degrees of freedom. Meanwhile, even small-area 
superconductor-insulator-superconductor (S-I-S) Joseph- 
son junctions, as solid-state mesoscopic systems, have 
many low-energy degrees of freedom inside amorphous 
insulating layer. Generally, those are phonons (i.e., vari- 
ables with oscillator-type equidistant energy spectrum) 
and the two-level (TL) systems (i.e., systems with non- 
equidistant energy spectrum so that only excitations 
from the ground state to the lowest excited state are in- 
volved at a given frequency). Previously it was shown 
that optical phonons inside intrinsic Josephson junctions 
in cuprate layered superconductors cause anomalies in 
the DC current-voltage characteristics. Namely, peaks 
in the tunneling current at the voltages corresponding 
to the phonon frequencies, V = fkJ p h/2e, were observed 
and mechanisms of their coupling with the junction phase 
difference were identified, see Refs. 1-4. 

While low-energy acoustic phonons are coupled weakly 
with the phase difference due to small statistical weight 
and negligible optical activity, the TL systems may be 
much more effective [5], A TL system may be, for ex- 
ample, an ion having two possible positions inside poten- 
tial wells with tunneling between them. These degrees 
of freedom interact with the junction phase difference if 
they modulate the Josephson energy and/or, if they are 
optically active. In the latter case they are coupled to the 
phase difference, 4>(t) (we measure the phase difference 
in units of magnetic flux) , via the electric field inside the 
junction, £ — <p/d, where d is the effective thickness of 
the junction. 

Here we propose a method to characterize the TL sys- 
tems inside Josephson junctions. Similar to phonons, the 
TL systems can cause anomalies in the DC I-V charac- 
teristics. But what is more, a TL system can precess 
at its Rabi frequency Or. when the Josephson oscillation 
frequency matches the level splitting, huj] = 2eV — HClo. 
This is similar to well known behavior of atom in the 
presence of a resonant electromagnetic wave, see Ref. 6. 



However, unlike the optical systems, where the Rabi os- 
cillation are observed as sidebands (f2o±flii) in the emis- 
sion, in our system, the Josephson current oscillates at 
the Rabi frequency itself. 

The transport and noise properties of the current and 
voltage biased Josephson junctions have been a subject 
of extensive studies. The techniques range from the semi- 
classical Langevin equation [7] which was generalized us- 
ing the quantum effective action [8] , to the infinite series 
summation [9, 10]. For V ^> I C R the Coulomb Blockade 
regime is realized, where the Cooper pairs tunnel inco- 
herently and give rise to the shot noise. We will consider 
this regime for a point-like Josephson junction (its size 
much smaller than Josephson length). We will also as- 
sume that 2eV ^> k-gT and that wj is higher than all 
the frequencies at which we want to measure the voltage 
noise. 

First we consider a junction with the Josephson energy 
modified by the interaction with TL system, described by 
the pseudospin-1/2 operator S 
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where Ej = <&qI c /2h, I c is the Josephson critical current 
without TL system, and $o = h/(2e) is the flux quan- 
tum. Coupling constants j = (j x , j y , j z ) characterize the 
modulation of the Josephson critical current by the TL 
system. We assume a small shunt resistance, R <C Rn, 
where Rn is the junction's normal state resistance. This 
ensures that voltage V on the junction can be small com- 
pared with the superconducting gap A in the running 
phase regime. The Hamiltonian of the voltage biased 
system reads [9, 11] 



H 
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where -ffR.(£) is the Hamiltonian of the resistor on which 
the phase £ drops, and q is the charge on the capacitor 
C, q being conjugate to (j). After the transformation <j> = 
4> — Vt and H = H — Vq we obtain 



H = 



2C 



H R {-4>) + Hj(4> + vt)-nn s z 



(3) 
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where q = q — CV. The Josephson current operator is 



The Josephson current operator transforms as 



I 3 =7 c (l+j-S)sin 



2ir(4> + Vt) 
$ 



(4) 



For frequencies such that hui <C 2eV and in the run- 
ning phase regime, V ^> I C R, the symmetrized phase 
autocorrelator S<f,(uj) is related to the symmetrized cur- 
rent autocorrelator, Si(uj), via 



"5,4,0 + 



i? 2 



Si, S$ 



HR 
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where = 1 + C 2 R 2 lo 2 . In Eq. (5) the equilibrium 

correlator S^o is due to the first two terms of Hamil- 
tonian (3), while the next (57) term accounts for the 
Josephson coupling. Eq. (5) can be obtained either from 
an exact relation between the phase and current Green's 
functions or from the quasiclassical Langevin equation. 
In typical experiments, RCoj <C 1 for all relevant fre- 
quencies uj, and, therefore, Y(u>) w 1. We, however, keep 
the factor Y(lo) to be able to discuss the regime of a 
junction shunted by a big capacitor. 

Next, we calculate the correlator of the Josephson cur- 
rent operators. We are particularly interested in a near- 
resonance situation, O Riwj. We will use the full quan- 
tum Hamiltonian (3) to calculate the spin's contribution 
to the correlator Si. As we can see from Eq. (3), the spin 
is subject to an ac driving at frequency u)j, "broadened" 
by the fluctuating phase <fi(t). Thus it is convenient to 
transform to the frame rotating with the angular velocity 
(2ir/$o)[V+(d/dt)(f)]. Formally this amounts to perform- 
ing canonical transformation H' = \JHU~ 1 + ihUU -1 , 
with 



U = exp 



^(4> + vt)S z 



(6) 



Without lost of generality we take j = (j±,0,j\\). The 
result is 

h'= |U h r (-4>) - h(n - uj 3 )s z - 

-gj(l+i||5,)cos 27r( '^ + ^ ) - m K S x , (7) 

5*0 

where Or = j^_I c /(Ae) is the Rabi frequency of the spin. 
The counter-rotating term (oc exp ± 4ir((j) + Vt) / $ ) can 
be shown to be not important. The resonance is reached 
when wj = Oo. Then the spin rotates around the x-axis 
(of the rotating frame) with the Rabi frequency Or , but 
its dynamics is affected by the noise due to the charge 
and phase fluctuations. 

The operator of the charge on the capacitor is trans- 
formed in the rotating frame as 



UqU~ 



2eS z 



(8) 



I'j = UIjU- 1 

= 4(1 + j|| S z ) sin 



2ir(j> + Vt) j ± I c 



Sit 



(9) 



which shows that dynamics of the spin translates into 
dynamics of the current. 

For |j| <C 1, the spin's contribution to the average 
current is negligible (unless spin's energy relaxation is 
much faster than that of the Cooper pairs). This, how- 
ever, is not always the case for the current near the 
Rabi frequency as described by the correlator Si (to) = 
(Ij(t)Ij(t'))u atwsi flu,- Indeed, as we shall see, the 
spin-dependent part of 7j gives rise to a singular (peaked) 
contribution to the correlator Si. We calculate Si(lu) in 
the rotating frame using Eqs. (7) and (9). The smooth 
part of Si(uj) is dominated by the shot noise of the 
Cooper pairs which tunnel incoherently The lowest or- 
der in RI c /[VY(u))] approximation gives good results for 
Sf ot if VY(uj) > I C R, uj < uj and R <C Rq = h/4e 2 . 
Then Sf ot {uj) = {I 2 /A)[P(lu + ujj) + P(uj - ujj)}, where 
P(t) = exp[J(t)], and J(t) = (27r/$ o ) 2 ((^(i)-0(O)) 0(0)) 
[12]. At high voltages 2eV S> k B T we obtain for sym- 
metrized correlator 



Sf ot (Q R ) « S} 



hot 



(0) « £p( Wj ) 



eRI 2 
VY(uj) 



(10) 



In addition, there exists a peak-shaped contribution 
to Si(u) near uj w Or. Using (sin27r(</> + Vt)/$o) = 
RI c /[2VY(ojj)] = I av we take into account the jy term 
by defining S* = j±S y + (I av / 'I c )j\\S z . Then 



S f«(t,t') = -?-{{S*(t)S*(t>)} + ) 



(11) 



To calculate this correlator we use Hamiltonian (7). Ex- 
actly at resonance the effective magnetic field M2r is di- 
rected along the x axis. Assuming the Rabi oscillations 
are under-damped (to be checked for self-consistency) we 
obtain 



gspin 



i 2 I 2 
16 



(w - Or) 2 + r 2 



+ (uj — » -uj) 



, (12) 



where j 2 s = j\ + (I av /I c ) 2 j 2 v 

To calculate the dephasing rate T2 we note that in (7) 
both the voltage noise SV = q/C and the shot noise (the 
j'll term) are coupled to S z , i.e., transversely to the Rabi 
field. Thus they contribute to T 2 through the longitudi- 
nal relaxation rate Ti, with T 2 — (1/2^. 

The voltage noise consists of two (uncorrelated to the 
lowest order) contributions: the equilibrium Johnson- 
Nyquist one and the one due to the shot noise of the 
Cooper pairs. Using 5V = (d/dt)<p (when no spin is 
present) and the expression for 5,4,0, we obtain for the 
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Johnson-Nyquist noise 



Sy (Or) = TTTT^ 7 COth ■ 



2k B T 



(13) 



The corresponding rate is I^ N = (e/h) 2 S^ (Or). 

The shot noise contributes twice: as part of the voltage 
noise SV and via the j|| term in (7). Using Eq. (5) for 
the SV part we obtain 



sf ot (n R ) 



R 2 



sf ot {n R ) . 



Y(n R ) 



(14) 



The corresponding rate is rf ot ' 5v = (e / 'Tif S*^ hot (fi R ) . 
Finally the j|| term contributes the rate 



rjj = (j||/4e) 2 S? hot (Q R ) 



(15) 



The last equation follows from ((sin tp(t) sin y(t'))) = 
({cos tp(t) cos <p(t'))), where = 2n((/>(t) + Vi)/$ - 

One can check that all three rates are to be added "in- 
coherently" , that is the noise cross-terms vanish. In par- 
ticular, the cross-term between the two shot noise contri- 
butions vanishes because ((sin (£>(£) cos <p(t'))) = 0. Thus, 
accounting for intrinsic TL dephasing Fq, 



-p -pJN i -pshot,(5V -pi i t 

1 2 — 1 2 + 1 2 +J- 2 + i 0- 



(16) 



For the "signal", i.e., the height of the voltage peak we 
obtain 



a1 r>2 j2 

s v eak (n R ) = J^.' 



16T 2 Y(n R ) ' 



(17) 



while the "noise", i.e., the smooth background is given 

by 



bg _ R 2 
bv " Y(n R ) 



R 



2k R T 



(18) 



Collecting all the terms we obtain the signal-to-noise ra- 
tio 



n = 



Y(n R )A, 



coth 



+ 



eR 2 I?B,, 



2k B T hfl R VY(uj) 



+ 



2R Q Y(n R )r 

tt-RQr 



(19) 



and R Q = h/(4e 2 ). 

For purely transverse coupling, ju = 0, the essential 
physics is the following: we illuminate the spin with the 
"magnetic" field 2ffi R cos 2ir(Vt + 0)/$ O - This field can 
be thought of as having a sharp peak (a line) near u) = u>j. 
The width of this Josephson line is given by the total 
voltage noise at zero frequency, hAv = ttSv(0)/Rq = 
(tt/R q )(S v n (lo = 0) + R 2 Sf ot (cj = 0)). This relation 
between the width of the Josephson line and the total 



voltage noise was obtained in Rcfs. [13, 14]. The Rabi os- 
cillation produced by this "line" are, in turn, also broad- 
ened by the same amount (in addition to the intrinsic 
broadening) T 2 = Alo/2 + r . Finally, the spin's (broad- 
ened) Rabi precession leads to broadened oscillations of 
the Josephson current and voltage at Or on top of the 
background of the JN and shot noise. 

It is also important to note that the Rabi oscillations 
of the pseudospin correspond to exactly one Cooper pair 
going back and forth across the junction. This can be 
seen from Eq. (8), or from the fact that for the Rabi os- 
cillations to occur exactly one "Josephson photon" with 
the energy Hloj must be absorbed and reemittcd by the 
spin, i.e., exactly one Cooper pair must go through the 
voltage drop V. 

To get a feeling for the relevant numbers we take the 
data obtained by Simmonds et al. [5], where the two 
lowest levels of a junction (phase qubit) in the super- 
conducting (phase-non-running) regime were driven res- 
onantly. The level splitting Woi(-T) was varied by the 
bias current / in the frequency interval 8.6-9.1 GHz. At 
some values of a>oi (-0 appreciable splittings (avoided level 
crossings) were observed. This was suggested to originate 
from TL systems with Oo ~ cjoi(I), and the interaction 
with the phase difference of the type (1). The splitting 
is caused by the j± term, while j\\ term is inessential as 
long as j\\Ej <C TiOo- Thus, the strongest impurity had 
j± w 6.5 • 1CP 5 , while we have an upper bound for the 
strength of the longitudinal coupling < 1CP 3 . This 
gives R w 2tt x 200 MHz (the spilitting of 25 MHz in [5] 
is due to the reduction factors corresponding to the zero- 
point motion of the phase degree of freedom in the poten- 
tial well). In Ref. [5] the critical current is I c w 10 /uA, 
the normal resistance of the junction is Rn ~ 30 O, while 
C ~ 1 pF. Using these parameters we estimate the signal 
strength and the signal-to-noise ratio 1Z. For the temper- 
ature we assume T w 10 mK, or k B T/h = 2tt x 200 MHz. 
The minimum voltage is given by I C R. We assume the 
shunting resistance of order R ~ 0.1O <C Rn- Shunts 
of this magnitude have been used in [15]. Hence, we 
have wj > (2e/h)I c R w 2tt x 0.5 GHz. From above wj 
is restricted by the gap which gives cjj < (2e/K)I c RN ~ 
2tt x 150 GHz. Thus we can take wj - 2tt x 10 GHz 
to be in resonance with the observed TLS. For j\\ = 
and assuming r = we obtain the signal-to-noise ratio 
1Z 0.25. For the maximally allowed = 10~ 3 the ratio 
1Z does not change considerably. 

For the integrated signal (signal amplitude) we obtain 

fo !^& 0ak l V2 w ?f RI < w 10- 2 nV. The Rabi line 

width is dominated by the Johnson-Nyquist noise, T 2 ~ 
2ir x 5 kHz. If Fo exceeds considerably this value the 
signal-to-noise ratio TZ will be reduced. 

We also note that for the above introduced parameters 
we have 1/{CR) ~ lO^s" 1 . Thus F(0 R ) w Y(wj) w 1. 
One has to increase C by at least three order of magni- 
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tude in order to start having Y(Or) > 1. 

For some (extreme) parameters, such that Y(Or) >• 1, 
we obtain 1Z » 1 which is in contrast with the limita- 
tion 1Z < 4 found for the measurements of the peak in the 
current noise at the frequency Oo [16-18] using a normal- 
state tunnel junction (broad band amplifier). In that case 
the voltage V ^> Kl^je (broad band) is applied. It in- 
coherently excites the TL system but also introduces the 
relaxation due to dissipation necessary for measurement 
procedure. The relaxation is determined by the noise at 
frequency flo and the signal is measured on the back- 
ground of the noise at the same frequency. As a result, 
1Z is a universal number. In the case considered here, 
spin is excited at (high) frequency Oo, but the signal is 
observed at low frequency, due to nonlinearity of coupled 
spin and Josephson junction system. It is worth mention- 
ing, however, that in the regime with large signal-to-noise 
1Z the single Cooper pair mainly charges and discharges 
the capacitor, barely going through the resistor; thus, the 
integated signal in this regime is reduced. 

Now let's consider the mechanism where spin couples 
to the junction via the electric field. The Hamiltonian is 



H 



2C 



+H K (Vt-cj>)-Ejcos 



2n(f> 
$ 



- TiQqS z 



C 



S x > 
(20) 

where Qtl is the effective charge of the TL system, given 
by Qtl = (Itl/L, where c?tl is the TL's dipole moment, 
while L is the junction's width. For simplicity we as- 
sumed purely transverse coupling. Remarkably, the split- 
ting observed in Ref. [5] could also be explained by (20) 
with Qtl ~ e. In this mechanism the spin is coupled to 
the variable q, which in zeroth order in tunneling (£j) 
does not oscillate and has only the Johnson-Nyquist noise 
spectrum. At V > RI C this variable acquires an oscillat- 
ing part due to the Josephson oscillations of the current. 
Thus the Rabi driving becomes possible. The width of 
the Rabi line is again determined by the full width of the 



Josephson line Ij +T 



nshot.5V 



From the integral (weight) 



of the Josephson peak in the S q (u>) correlator we obtain 
the Rabi frequency Or = RI C Qti,/ {2fiyjY{ujj)). Then, 
analysis similar to the one presented above again gives 
Eq. (19) for the signal-to-noise ratio (with A\\ — B\\ — 1 
as we assumed purely transverse coupling). For Qtl ~ e, 
we obtain a similar to the previous mechanism Or and 
a similar value of 1Z. For the integrated signal we obtain 

1/2 



^peak 



QtlA 2 Jc 
eR Q 



10" 2 nV. Note that the 



f dul q 

Rabi frequency Or and the integrated signal depend dif- 
ferently on R in two coupling mechanisms. This may al- 
low to distinguish between the two, while they are undis- 
tinguishable in measurements of type [5] . 

In this letter we discussed what happens when the 
Josephson oscillations are in resonance with one TLS. Let 
us mention another interesting possibility to manipulate 
the system. By changing the applied voltage slowly, one 
can create the regime of the "adiabatic passage" when 



u)j(t) passes slowly via Oo and exactly one additional 
Cooper pair is transferred through the junction. Varying 
u)j(t) in a wide enough interval one can "touch" many 
TL systems and create a measureable additional current. 

In conclusion, we propose that the measurements of 
the low frequency voltage noise in a Josephson junction 
in the dissipative (running phase) regime may be used 
to characterize the TL systems inside the junctions, i.e., 
energy splitting O , coupling strength j± from the Rabi 
frequency, and intrinsic dephasing rate r from the height 
of the voltage peak, Eq. (17). We predict a peak at the 
Rabi frequency when a TL system is resonantly driven by 
the Josephson oscillations, wj = Oo, with the Rabi fre- 
quency proportional to the interaction strength between 
the TL system and the Josephson phase. The peak inten- 
sity (signal-to-noise ratio) can be controlled by the shunt 
resistor and capacitor. 
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